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Fermi-edge singularity in the vicinity of the resonant scattering condition 

V. V. Mkhitaryan and M. E. Raikh 
Department of Physics, University of Utah, Salt Lake City, UT 84112, USA 

Fermi-edge absorption theory predicting the spectrum, A(lij) oc ij~'^°oi'^+ o/'^ ^ relies on the as- 
sumption that scattering phase, 5o, is frequency- independent. Dependence of So on u becomes 
crucial near the resonant condition, where the phase changes abruptly by tt. In this limit, due to 
finite time spent by electron on a resonant level, the scattering is dynamic. We incorporate this 
time delay into the theory, solve the Dyson equation with a modified kernel and find that, near the 
resonance, A{u>) behaves as ijj~ ' \\nuj\. Resonant scattering off the core hole takes place in ID and 
2D in the presence of an empty subband above the Fermi level; then attraction to hole splits off a 
resonant level from the bottom of the empty subband. Fermi-edge absorption in the regime when 
resonant level transforms into a Kondo peak is discussed. 

PACS numbers: 74.25. Gz, 74.50. +r,73.40.Gk 
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Introduction. Many-body character of absorption from 
the localized level into continuum in the presence of the 
Fermi sea [l| manifests itself via two mechanisms: scat- 
tering of the excited electron from the hole left behind 
and adjustment of the Fermi sea to the abrupt switch-on 
of the hole potential. Correspondingly, near the thresh- 
old, the spectrum, A(uj) = I{uj)P{uj), is a product [2| of 
two functions 
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(1) 



where D is the bandwidth; ^o is the phase shift for s- scat- 
tering. Phase ^0 is related to the hole potential strength 
V as tan So = tti/qF, where vq is the density of states. 
The function /(w) is the single-electron-line contribution, 
which describes the first mechanism, while P{i-u) accounts 
for the shake-up effects. 

Derivation of Eq. ^ by Nozieres and De Dominicis 
(ND) in Ref. [4I strongly relies on the assumption that Sq 





FIG. 1: (Color online) (a): In the vicinity of a topological 
transition, the core hole potential creates bound state (reso- 
nant level) with energy, Eex, measured from the bottom of 
the upper subband; its lifetime, F, is due to a coupling to 
the lower subband. For a photoelectron with energy, uj <^T 
above the Fermi level, dynamical character of scattering from 
resonant level is crucial, (b): Contribution to the absorption 
amplitude which is strongly affected by resonant character of 
scattering: two scattering acts involved, (2 + 3) and (4 + 5) , 
have opposite signs of amplitudes, Eq. (llUp . 



is frequency-independent. This assumption implies that 
the scattering is instantaneous. In the present paper we 
focus on the situation when this assumption is violated. 
This situation realizes when the Fermi level is close to 
an empty upper subband, as shown in Fig. [1] In such 
an arrangement it is important that even a deep core- 
hole creates a localized level a distance E^x below the 
bottom of the upper subband, see Fig. [T] Formation of a 
localized level happens in low dimensions (2D and ID). 
Experimental and theoretical studies of the arrangement 
in which Fermi level is near the bottom of empty subband 
have been reported in [3|46| and |7H13j. respectively. 

Due to degeneracy with continuum of lower band, the 
level Eex acquires a finite width, F. Then the process un- 
derlying the Fermi-edge singularity is not scattering from 
bare core-hole potential but rather resonant scattering 
on the quasi-local level. Correspondingly, the scattering 
pase is given by 

r 

S(uj) = arctan -, (2) 

where A is the energy distance from localized level to the 
Fermi level, Ep. In the limit A — ;■ 0, w-dependence of S is 
strong. Indication that in this limit the ND theory is in- 
applicable follows from the fact that at w = A, the phase 
abruptly changes by tt. This, in turn, leads to a paradox- 
ical conclusion that, as Fermi level is swept through the 
bound state, /(w) changes abruptly from I{uj) oc 1/w to 
I{uj) oc uj. By contrast, the shake-up term P{uj) oc w^/* 
does not experience a jump. 

As a resolution of this paradox, in the present paper 
we show that abrupt change of /(w) does take place in 
the frequency interval, lo '^ |A|, which gets progressively 
narrow as A goes to 0. We also show that outside this 
interval, where the ND theory does not apply, a new, dy- 
namical resonant-scattering regime governs the absorp- 
tion [14|. We show that in this regime, instead of Eq. 
((H), I{uj) is given by 
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FIG. 2: (Color online) Absorption spectrum, A{ij), in the 
situation when the Fermi level is close to resonant level 
(|A| <^ r) is shown schematically. For oj < |A| the u- 
dependence, Eq. (|29l) . is governed by ND theory. In the 
"dynamic scattering" domain |Aj < oj < F the spectrum is 
A{ui) = (F/lj)^/* ln{Ee^/uj) ln(Ee^/r) for any sign of A. 



independently of the sign of A. Our findings are illus- 
trated in Fig. [21 

Dynamic-scattering regime. Within the ND theory, the 
transient Green function, (p{T, r'ji, t'), satisfies the Dyson 
equation 

^{r, T'\t, t') = G{t - r') + [dTK{T, f )^(f , r'|i, t'), (4) 



where G{t) = —v^jt is the bare electron Green function, 
and the kernel is defined as K(t., f) = VG(t — f). The 
single-electron-line contribution to Aiyi) is expressed in 
terms of the transient Green function as 

1 fO 

/(cj) = Im/ dte"*"V(0",i+|i,0). (5) 

Delicate character of Fermi-edge absorption can be in- 
ferred from the structure of n-th term. 



F"/ dtx' 



dt^ G{-ti)G{ti - i2) • • • G(t„ - i), (6) 



of perturbative expansion of iy9(0~, t+|t, 0) in powers of 
V, which follows from Eq. (|3]). It appears that contribu- 
tions of different orderings of times ti , Fig. [5^, at which 
electron is scattered, cancel each other up to 1/nl. 

Frequency dependence of S leads to a time delay, r, in 
the scattering processes, as illustrated in Fig. [3}d. Each 
delay ranges from E~^ to |A|^^. As a result, in the 
limit A — )■ 0, the above cancellation is completely de- 
stroyed. On the quantitative level, dynamical character 
of the resonant scattering alters the kernel of Eq. ^ . 
Resonant- scattering kernel. In the presence of the empty 
upper subband, the constant interaction V for electrons 
of the Fermi sea acquires a frequency-dependence of the 
form 



^(w) = Vr. 
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FIG. 3: (Color online) (a): Typical contribution, Eq. ((6|, 
to the transient Green function in ND theory. Instantaneous 
scatterings take place at time moments, fi, t2,...,t„. (b): In 
the presence of a resonant level, each scattering act involves 
a "visit" to the resonant level, and takes a finite time, r. 
Blowup illustrates a ladder of scatterings off the core-hole in 
the upper subband, 2, which lead to resonant scattering in 
the lower subband, 1. 



Here V12 and V22 are inter- and intra-subband matrix 
elements of the hole potential, respectively; G2 is the 
Green function of the upper subband, 

1 



G2(^) 



E 



uj - e2q + iri 



(8) 



and e2q — h^q^/2m is the spectrum near the bottom. 
The form Eq. ([7]) is a result of summation of ladder-type 
diagrams, see a blowup in Fig. |31 Bound state emerges 
as a pole in Eq. d?]). 



E, 



2TT^ni 
~l9~ 



T/22P (ID), E,^ = De^^\^ (2D). (9) 



Expanding near the pole, we simplify V to the form 

V{uj) = (—] 1--, (10) 

where width, F, is given by 

r = ^-^\v,2? (ID), r = ^-^^^^\v,2? (2D). 

(11) 

Effective interaction Eq. (jlOp corresponds to repulsion 
for w > A and to attraction for w < A. In the resonant- 
scattering regime, the kernel in Eq. ^ instead of the 
simple product VG{t — f) becomes a convolution 



Kres{T,T')^ dx* G{t - T 



duj 
2^' 



lUJiT —r 



V{u^) 



(12) 
Integration over uj can be carried out explicitly, yielding 



Kres{T,T') = i- 



dT* 



(ll_llO^-jA(r*-r') 



(13) 



For most interesting case, A = 0, the expression for the 
kernel assumes the form 



Kres{T,T') = -i-lll 



(14) 



We see that In |t — r'| in Kres emerges in the place of 
l/(r — r') in the ND kernel. This is a consequence of 
opposite signs of interaction with resonant level for elec- 
trons and for holes. 

Derivation of Eq. ^. Substituting Eq. p^ into Eq. 
(HI and performing the rescaling 



^(r,r'|i,0) 



t 
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(15) 



we arrive to the following dimensionless equation 



(t>{x,y) 



1 



x~y 



ij / dz In 



{z,y), (16) 



where 7 = —t. Note that the kernel Eq. p6)) is much less 
singular than the ND kernel. As a result, the power-law 
transient factor in the ND solution, [{l—x)y/{l—y)xY^^, 
which is singular for a; — )■ 0, y — )■ 1, does not emerge. This 
allows us to search for solution of Eq. (IT51) as a linear 
combination, 



[x,y) 






Cn{y,l)Un{x), 



(17) 



of the eigenfunctions, Un{x), of the Hermitian integral 
operator 



R{Un{x}} = dz\n\x - z\Un{z) ^ XnUn{x), (18) 

Jo 



where A„ are the eigenvalues. The kernel in Eq. (|T8|) 
is symmetric with respect to x = 1/2. Correspondingly, 
Un{x), n — 0, 1, 2, • • • are even (for n even) and odd (for 
n odd) functions of (x — i). Substitution of Eq. (|17p into 
both sides of Eq. ([16]) yields the following expression for 
coefficients of expansion. 



liy,!) = 



am{y) +b.md{y,-f) 

1 + i-fXm 



(19) 



where 



ani{y) = / dz-^ = Xmu'^iy), (20) 

Jo y~ z 

bra= dzUm{z)\nZ^XjnUra{0). (21) 

Jo 

In the second identities of Eqs. (|20p . (I^TI) we used the 
properties of operator R. 

The kernel, Inx — In jx — y|, of Eq. (|16|) . in addition 
to the difference term, contains a y- independent term. 
As a result, bm in the expression for Cmiy, 7) enters with 
coefficient 



d{y,7) =«7 / dz(l>iz,y), 
Jo 



(22) 



which is the same for all m, so that Eq. (122]) can be 
viewed as a self-consistency condition. 



Our key observation which will be justified later is 
that, the relevant values of m in the expansion Eq. (J17p 
are ^ 1. For such m, it can be shown that eigenval- 
ues behave as A^, ~ —1/m. Concerning the eigenfunc- 
tions Um{x), they assume constant value near the bound- 
aries, Um{0) — ±Wm(l), while the derivatives, u'^{x) di- 
verge logarithmically at the boundaries. Outside of small 
(^ 1/m) intervals from the boundaries, they behave as 
sin(27ra;/Am) and cos(27rx/Am). AH these large- m prop- 
erties can be established upon integrating by parts in Eq. 
(1181) . Then the remaining integral will be determined by 
a narrow domain, |a; — z| ^ l/m. The same simplifica- 
tion allows us to obtain a concise expression for d{y, 7), 
namely, 



diy,"/) ^ -i7ln(l-y). 



(23) 



and subsequently the solution for 0(a;, y) in the form 

[u„(l) + ijXnUn{Q)]Un{x) 



{x,y) ~ -ln(l -y) 



1 + i"fX„ 



(24) 
In calculating /(w) from Eq. ((M)) the combination, 
^„ M„(l)u„(x), appears. This combination is equal to 
(5(1 — x) and does not contribute to absorption, which 
is determined by x — > and ?/ — > 1. In ND problem, 
divergences in the time domain are cut off by the inverse 
bandwidth, iD~^. In our case, the logarithmic diver- 
gence in Eq. (f24|) is terminated at (1 — y) — i{Eext)~^ ■ 
Correspondingly, the minimal x in Eq. (j24l) is, in fact, 
Xmin — i{Eext)~^ ■ Taking this into account and using 
Eqs. (O and ([T5|, we express /(w) in the integral form 



I{ijj) ~ Im 



dte''^^ S{t)\n{iE^xt), 



where the sum, 



Sit) 



3dd m 



XiriXt 



- ITT 



(25) 



(26) 



is performed over only odd m, for which Um{l) = 
— Um(0). Summation over m requires large- r/i values of 
Um(0). In fact, these values saturate with increasing m. 
This follows from the identity 



E 

odd m 



1 ( X 

X,nU,n{Q)u,n{x) = - In 

2 V 1 — a; 



(27) 



The left-hand side has a logarithmic divergence at a; — ;> 0. 
In the right-hand side, with Am ~ —l/m, logarithmic 
divergence of the sum is achieved with Umix) « Mm(0) ~ 
1 for m < x^^. Eq. ^7\ is a direct consequence of the 
identity, Ina; — J2m'^rn'U'm{0)umix). The sum over odd 
m appears upon writing this identity for x and (1 — x) 
and taking their difference. Comparing Eqs. (j26p and 
(P7)) we conclude that 



Sit) ~- + Tr- in(*rte, 

4 ZTTl 



4 + ii"(T^M^«> 



Imaginary part of S{t) is determined by the upper cut- 
off, while ReS'(i) comes from m ~ Ft; relevant values 
of m are large, as we assumed above. Factor Tt in the 
argument of the logarithm comes from the terms in sum 
Eq. (HSl) with m < Tt. 

With S{t) being time-independent for t S> F^^, inte- 
gration in Eq. (|25p recovers our result Eq. ([3]) . Physical 
meaning of the enhancement factor, F/w, in Eq. ([3]) is 
the number of times electron virtually visits the resonant 
level and returns back to the lower subband during the 
time l/uj. 

Concluding remarks, (i) Our main result, Eq. (|3]), was 
obtained by setting A = in the kernel Eq. ([T3|. 
Finite-A correction to A{uj) is small as (A/w)^. Let us 
briefly discuss the opposite limit, w ^ |A| ^ F. In this 
limit, the characteristic time t in Fig. [3] is ;» l^l~^- 
Then, our above consideration for A = applies only 
in the time domain (— |A|~^,0). In the remaining do- 
main (i, — |A|~^) characteristic time intervals ti — t^+i in 
Fig. |3]are much bigger than |A|~-'^. This makes scatter- 
ing " instantaneous" , and ND theory applicable in the do- 
main (t, — |A|~^); corresponding time-independent scat- 
tering amplitude and effective bandwidth are T/{'!ti>qA) 
and A, respectively. Separation of time scales leads 
to the absorption coefficient in the form of a product 
of (a;/|A|)-2arctan(r/A)/,r ^ (a;/| AD-'^'S"^ (from long 
times), and F/|A| (from short times, \t,\ < \A\-^). With 
shake-up included, the overall result for absorption coef- 
ficient, 



A(^) 



F 

oc 

A>0 U! 



F /W\ 4 . , , WF /W\ 4 , , 



is illustrated in Fig. [2j At w '- A, Eq. ([29l) matches Eq. 
([3]) with logarithmic accuracy. 

(ii) Ellipses in Fig. [3] can overlap; this corresponds to 
double-occupancy of the resonant level, which our the- 
ory does not prohibit. This implies that the Hubbard 
repulsion, {/, is smaller than F. Since U is inversely pro- 
portional to the radius of the level wave function, this 
condition can be met. 

Even more interesting is the situation when resonant 
level, which is split off the upper band by the core-hole, 
is well below the Fermi level (|A| ^ F), and U is not 
small. Then the level gets occupied, and after time 
~ T^^ ~ exp[7rA/2F] a Kondo peak of width ~ Tr- is 
formed at the Fermi level ^7,] . Photoexcited electron will 
experience dynamical resonant scattering from this peak. 
One can argue that our result for dynamical resonant 
scattering regime will apply in this situation, with re- 
placement, (F/cj) — ^ (Tk/i-^). However, we cannot make 
this analogy rigorous because the Kondo scattering of the 
state u! leads to finite lifetime ~ Tk/oj'^. 

Another interesting possibility is then the core-hole 
possesses spin [Ig. Then a Kondo peak at the Fermi 
level will form even without the upper subband, so that 



the scattering off the core-hole will become a Kondo- 
scattering with w-dependent scattering phase. 

{in) We stated above that, unlike the single-electron- 
line contribution, I{oj), the shake-up factor, P{oj) = 
(^/P)arctan=(r/A)/7r=^ « (^/r)^/", is uot Sensitive to the 
sign of the detuning, A. To prove this statement rig- 
orously, one has to extend evaluation of the many-body 
overlap integral Ref. [19[ to the case of resonant scat- 
tering with frequency-dependent scattering amplitude 
V{uj). The modified overlap integral reads 



1 



(27r)2 



dX / dQ- 



v{n)Gin) 
'i~xv{n)G{n) 



dn 



,\n[i-xv{n')G{n')] 
{fi' ~n + M/)2 



(30) 



For V{n) = const(ri), the integral Eq. (130]) comes from 
small arguments, namely J7 — t- 0+ and 57' — > 0~. It di- 
verges logarithmically with coefhcient —{Sq/tt)'^, where 
So = — Arg[l — VGr{0)], and Gr is the retarded Green 
function. Using the analytical properties of resonant 
scattering amplitude, Eq. ([TU|. one can check that 
the above logarithmic divergence persists when V is 
frequency-dependent. Independence of the sign of A fol- 
lows from the fact that coefficient in front of logarithmi- 
cally divergent part is ;^ (5a(0)(5_a(0), where (5a (w) = 
tan~^[F/(a; — A)] is the resonant-scattering phase. 

(iv) We neglected a contribution to absorption from a 
direct transition core- level — >■ i^ex- Since the width of 
Eex is ~ F, this contribution is non-singular at small uj. 

(v) Fermi edge physics manifests itself also in resonant 
tunneling via an impurity J20j . However, the regime A <C 
F, considered in the present paper, can not be realized in 
this setting since at long times impurity gets reoccupied. 
Acknowledgments. We are grateful to T. V. Shahbazyan 
for very helpful discusssions. 



[1] 

[2] 

[3] 



[4] 
[5] 

[6] 

[7] 
[8] 



G. D. Mahan, Phys. Rev. 163, 612 (1967). 

P. Nozieres and C. T. De Dominicis, Phys. Rev. 178, 

1097 (1969). 

M. S. Skolnick, D. M. Whittaker, P. E. Simmonds, T. 

A. Fisher, M. K. Saker, J. M. Rorison, R. S. Smith, P. 

B. Kirby, and C. R. H. White, Phys. Rev. B 43, 7354 
(1991). 

W. Chen, M. Fritze, W. Walecki, A. V. Nurmikko, D. 

Ackley, J. M. Hong, and L. L. Chang, Phys. Rev. B 45, 

8464 (1992). 

T. Melin and F. Laruelle, Phys. Rev. Lett. 85, 852 (2000); 

T. Mehn and F. Laruelle, Phys. Rev. B 65, 195303 

(2002). 

K. Behnia, L. Balicas, and Ya. Kopelevich, Science 317, 

1729 (2007). 

L. Balents, Phys. Rev. B 61, 4429 (2000). 

M. Sitte, A. Rosch, J. S. Meyer, K. A. Matveev, and M. 

Garst, Phys. Rev. Lett. 102, 176404 (2009). 



[9] O. A. Starykh, D. L. Maslov, W. Hausler, and L. I. Glaz- 

man, Lect. Notes Phys. 544, 37 (2000). 
[10] N. N. Ablyazov, M. Yu. Kuchiev, and M. E. Raikh, Phys. 

Rev. B 44, 8802 (1991). 
[11] J. Alicea and L. Balents, Phys. Rev. B 79, 241101 (2009). 
[12] J. F. Mueller, Phys. Rev. B 42, 11189 (1990). 
[13] T. V. Shahbazyan, N. Primozich, I. E. Perakis, and D. S. 

Chemla, Phys. Rev. Lett. 84, 2006 (2000). 
[14] Dynamic and non-equilibrium aspects of conventional 

Fermi-edge singularity have been considered in Refs. [1^ 

and flo], respectively. 
[15] N. Primozich, T. V. Shahbazyan, I. E. Perakis, and D. 

S. Chemla, Phys. Rev. B 61, 2041 (2000); D. Porras, 

J. Fernandez-Rossier, and C. Tejedor, Phys. Rev. B 63, 

245321 (2001); N. dAmbrumenil and B. Muzykantskii, 

Phys. Rev. B 71, 045326 (2005). 



[16] T. K. Ng, Phys. Rev. B 54, 5814 (1996); M. Combescot 
and B. Roulet, Phys. Rev. B 61, 7609 (2000); B. Brau- 
necker, Phys. Rev. B 68, 153104 (2003); B. Muzykantskii, 
N. dAmbrumenil, and B. Braunecker, Phys. Rev. Lett. 
91, 266602 (2003); D. A. Abanin and L. S. Levitov, Phys. 
Rev. Lett. 94, 186803 (2005). 

[17] P. Nordlander, M. Pustilnik, Y. Meir, N. S. Wingreen, 
and D. C. Langreth, Phys. Rev. Lett. 83, 808 (1999). 

[18] V. V. Mkhitaryan, E. G. Mishchenko, M. E. Raikh, and 
L. L Glazman, Phys. Rev. B 80, 205416 (2009). 

[19] D. R. Hamann, Phys. Rev. Lett. 26, 1030 (1971). 

[20] K. A. Matveev and A. L Larkin, Phys. Rev. B 46, 15337 
(1992); D. A. Abanin and L. S. Levitov Phys. Rev. Lett. 
93, 126802 (2004). 



